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Abstract 

This work examined the wake aerodynamics of a sin- 
gle helicopter rotor blade with several tip shapes operat- 
ing on a hover test stand. Velocity field measurements 
were conducted using three-component laser Doppler 
velocimetry (LDV). The objective of these measure- 
ments was to document the vortex velocity profiles and 
then extract the core properties, such as the core radius, 
peak swirl velocity, and axial velocity. The measured 
test cases covered a wide range of wake-ages and several 
tip shapes, including rectangular, tapered, swept, and a 
subwing tip. One of the primary differences shown by 
the change in tip shape was the wake geometry. The 
effect of blade taper reduced the initial peak swirl veloc- 
ity by a significant fraction. It appears that this is ac- 
complished by decreasing the vortex strength for a given 
blade loading. The subwing measurements showed that 
the interaction and merging of the subwing and primary 
vortices created a less coherent vortical structure. A 
source of vortex core instability is shown to be the ra- 
tio of the peak swirl velocity to the axial velocity deficit. 
The results show that if there is a turbulence producing 
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region of the vortex structure, it will be outside of the 
core boundary. The LDV measurements were supported 
by laser light-sheet flow visualization. The results pro- 
vide several benchmark test cases for future validation of 
theoretical vortex models, numerical free-wake models, 
and computational fluid dynamics results. 


Nomenclature 

A rotor disk area = nR 2 , m 2 

BEMT blade element momentum theory 
Cj rotor thrust coefficient = T /p A ( Q.R ) 2 

Ce Ludwieg vortex stability swirl coefficient 

C z Ludwieg vortex stability axial coefficient 

c rotor blade chord, m 

/ Gudermannian function 

k[ ratio =T V /Tb 

kn ratio =(r t ? /QRc) / (CV /<?) 

/ mixing length 

Nb number of blades 

r radial distance from vortex core center, m 

r c vortex viscous core radius, m 

R rotor blade radius, m 

Re Reynolds number 

Re „ vortex Reynolds number = F v /v 

Ri Richardson number 

Ri re-laminarization parameter 

Ric re-laminarization parameter at r ~r c 

S Leibovich-Stewartson vortex stability param. 



S p Bradshaw’s shape parameter 

t 0 effective time offset based on r c at t = 0, sec 

t time, sec 

T rotor thrust, N 

TPP tip path plane 

V ax axial velocity deficit in the vortex core, m/s 

Vq vortex tangential (swirl) velocity, m/s 

y dimensional radial distance along blade, m 

z vertical distance from TPP, m 

a Iverson’s empirical constant = 0.01854 

tj similarity variable = r 2 / (4a 2 T v x) 

y vortex circulation profile, T /T v 

y tr transitional vortex circulation profile 

T local vortex circulation, m 2 /s 

T b peak bound circulation, m 2 /s 

T v total vortex circulation, m 2 /s 

F c vortex core circulation = 27 tV^ m r C9 m 2 /s 

X mean inflow velocity ratio through the TPP 

v kinematic viscosity, m 2 /s 

vj effective turbulent viscosity, m 2 /s 

Q. angular velocity of rotor blade, rad/s 

p air density, kg/m 3 

a rotor solidity = Nbc/%R 

x time variable = t 4- to 

£ wake age of vortex core from TE, deg 


Introduction 

Many experiments have been attempted to examine the 
formation of lift-generated tip vortices. Two of the 
most detailed surveys are given by Rossow (Ref. 1) 
and Spalart (Ref. 2). Some of the first phase-resolved 
measurements of helicopter rotor flow fields were at- 
tempted using hot-wire anemometry (e.g., Refs. 3-8). 
Modern optical flow diagnostic techniques, such as 
laser Doppler velocimetry (LDV) and particle image ve- 
locimetry (PIV), can alleviate many of the limitations of 
hot-wire anemometry. An early application of a one- 
component LDV system to a rotor flow field was by 
Landgrebe & Johnson (Ref. 9), who provided a demon- 
stration of its potential capabilities. Another application 
was that of Lorber et al. (Ref. 10) using 2D LDV to mea- 
sure the circulation field around the blade of a Mach- 
scaled UH-60A with surface-mounted pressure trans- 
ducers. 

Over the last thirty years, LDV has gradually matured 
to allow simultaneous three-component flow field mea- 
surements with excellent spatial and temporal resolution 
(e.g., Refs. 1 1-15). Recently, PIV has become a rapidly 
maturing and competitive technique, as shown by the 
measurements of Heinecketal. (Ref. 16) and McAlister 
et al. (Ref. 17). A summary of the relative capabilities 


and limitations of the LDV and PIV techniques, as they 
apply to the measurement of rotor vortex flows, is given 
by Martin et al. (Ref. 18). The current work builds on 
the measurements reported by Leishman and co-workers 
(Refs. 15, 19-21). 

Results of detailed flow visualization (Refs. 22 
and 23) have documented the viscous flow field sur- 
rounding the rotor blade tip. The physical size of these 
vortical flow structures is only a fraction of the blade 
thickness, and so a fine level of spatial resolution is re- 
quired to capture and resolve details of the vortex core. 
A fine level of temporal resolution (i.e., phase resolu- 
tion) is required by the dynamic nature of the interaction 
and convection of the vortex filaments. These resolution 
requirements are critical when comparing different tip 
shapes. 

The present work involves three-dimensional veloc- 
ity measurements acquired in the flow field surrounding 
the tip of a one-bladed rotor in hover. An analysis of 
the measurement uncertainties involved in making these 
types of measurements using laser Doppler velocimetry 
(LDV) is described in Ref. 24. As a result of this initial 
study, the LDV technique was significantly refined and 
enabled the measurements presented here. Examples of 
the preliminary and refined techniques are given for sev- 
eral of the test cases. The results are from Refs. 25 and 
26, and are compared to the baseline test case presented 
in Ref. 23. 

Phase-resolved LDV measurements documenting the 
formation, roll-up, and evolution of the blade tip vor- 
tex were acquired by making a radial traverse across the 
vortex core at various planes below the rotor. This al- 
lowed the core properties to be studied as a function of 
wake age. The measurements for three other tip shapes 
were compared with the baseline blade. The planforms 
included a 20° swept tip, a tapered tip with a 0.4 taper 
ratio, and a sub wing tip. The vortex properties estimated 
from the velocity profiles were compared to the baseline 
measurements. 

Physical Description of a Rotor 
Wake 

The combined wake of a hovering rotor can be visual- 
ized as a system of interacting helical surfaces of vor- 
ticity (vortex sheets) and discrete tip vortices imposed 
within a streamtube of actuated fluid. It is the rotational 
motion of the rotor blade about the shaft that causes the 
spanwise lift distribution on the blade to be biased out- 
board, toward the tip. A representative lift distribution 
is shown in Fig. 1. As a result, a strong static pressure 
gradient forms around the tip region, and the outboard 
portion of the wake surface experiences a roll-up into an 




Vortex core 


Figure 1 : Schematic showing the bound circulation dis- 
tribution and the corresponding trailing wake on a rotor 
blade. 

intense region of vorticity in the form of a trailing tip 
vortex. The near-field tip vortex aerodynamics involve 
the confluence of the trailing edge turbulent wake and 
any edge- separated flow into a viscous vortex core. As 
will be shown, this process is a function of tip shape. 

Figure 1 shows that a fixed azimuthal distance from 
the blade trailing edge is defined as the wake age, 
This figure also shows the qualitative relationship be- 
tween the bound circulation distribution on the blade 
and the trailing wake structure. For a multi-bladed ro- 
tor in hover, the details of the trailing wakes are shown 
by Fig. 2 in an axial/radial plane at a fixed value of wake 
age, It is evident that the wake generated by a heli- 
copter rotor is characterized by highly dynamic, inter- 
acting vortical flow structures (e.g., Ref. 27). 

In contrast to a fixed wing aircraft, the rotor tip vor- 
tices remain close to the blades. This is especially true 
under certain flight conditions such as low-speed de- 
scending flight, and during transition to forward flight. 
Not only does the wake interact with the main rotor 
blades, but the filaments may intersect the tail rotor, 
the vertical and horizontal stabilizers, and the airframe. 
Even in hover, the vortices pass within a fraction of a 
chord length below the following blades. This small 
blade/vortex miss distance results in high induced veloc- 
ities including both upwash and downwash components, 
as well as radial flow. The velocity gradients induced by 
an individual blade’s trailing vortex has been shown to 
cause stall on the tip of the following blade - see Clark 
and Leiper (Ref. 28). The magnitude and direction of 
these induced velocities depends on the vortex strength, 
the distance from the blade, and the orientation of the 



Figure 2: Schematic showing a cross-section of the wake 
of a two bladed rotor operating in hover. 

filament with respect to the blade. A primary example 
of this is shown by the measurements of Berenger et al. 
(Ref. 29). 

The proximity of the wake to the rotor in hover creates 
a condition where small perturbations of the rotor inflow 
can cause the vortices to remain in the plane of the rotor 
during descending flight, in transition to forward flight, 
and during transient gusts. As will be shown by this fun- 
damental set of experiments, the effects of both swept 
and tapered tips are to keep the vortices closer to the 
TPP than for a rectangular tip in hover. 

As a direct result of the vortices remaining near the 
TPP, substantial unsteady variations in the local in- 
duced velocity field are then encountered by the rotat- 
ing blades. Blade- vortex interaction (BVI) manifests as 
a rapidly fluctuating angle of attack that results in an un- 
steady local blade loading. One key to predicting the 
magnitude of these angle of attack variations is to de- 
velop a better understanding of the vortex core structure 
(see Ref. 26). It is the coupling of the wake induced 
velocities to the rotor blade dynamics that makes un- 
derstanding the wake physics important for steady flight 
predictions, and critical for maneuvering flight predic- 
tions and/or simulations. 

Experimental Apparatus 

Rotor Test Stand 

The planforms tested in these experiments included a 
20° swept tip (sheared), tapered tip with a 0.4 taper ra- 
tio, and a subwing tip. The blade parameters and test 
conditions are given in the Appendix. The rotors were 
tested in the hovering state in a specially designed flow 
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conditioned test cell. The volume of the test cell was 
approximately 12,800 ft 3 (362 m 3 ) and was surrounded 
by honeycomb screens. This cell was located inside a 
large 500,000 ft 3 (14,000 m 3 ) high-bay laboratory. The 
rotor disk area was 5.6 ft 2 (0.52 m 2 ), and the wake was 
allowed to exhaust approximately 18 rotor radii before 
encountering flow diverters. Aperiodicity levels in the 
rotor wake when tested in this cell have been measured 
and used to correct the velocity field measurements (see 
Ref. 26). 

LDV System 

A fiber-optic based LDV system was used to make three- 
component velocity measurements. A beam splitter sep- 
arated a single multi-line Argon-Ion laser beam into 
three pairs of beams (green, blue and violet), each of 
which measured a single component of velocity (see 
Fig. 3). A Bragg cell set to a frequency shift of 40 MHz 
produced the second shifted beam of each beam pair. 

The laser beams were passed to the transmitting optics 
by a set of fiber-optic couplers with single mode polar- 
ization preserving fiber optic cables. The transmitting 
optics were located adjacent to the rotor, and consisted 
of a pair of fiber optic probes with integral receiving op- 
tics, one probe for the green and blue pairs, and the other 
probe for the violet pair. Beam expanders with focusing 
lenses of 750 mm and beam spacers were used to de- 
crease the effective size of the measurement volume. 

To further reduce the effective size of the probe vol- 
ume visible to the receiving optics, the off-axis backscat- 
ter technique was used, as described in Martin et al. 
(Ref. 24) and Barrett & Swales (Ref. 32). This tech- 
nique spatially filters the effective length of the LDV 
probe volume on all three channels. Spatial coincidence 
of the three probe volumes (six beams) and two receiv- 
ing fibers was ensured to within a 25 jam radius using 
an alignment technique based on a laser beam profiler, 
as described in Ref. 24. Alignment is critical for 3- 
component systems because it is geometric coincidence 
that determines the spatial resolution of the LDV probe 
volume. In the present case, the resulting LDV probe 
volume was measured to be an ellipsoid of dimensions 
80 jum by 150 ju m, which was about 3% of the maximum 
blade thickness (or about 0.5% of the blade chord). 

The flow into the rotor was seeded with a thermally 
produced mineral oil fog. A mineral oil based fluid is 
broken down into a fine mist by adding nitrogen. The 
mist is then forced into a pressurized heater block and 
heated to its boiling point where it becomes vaporized. 
As the vapor escapes from the heat exchanger nozzle it 
is mixed with ambient air, rapidly cools, and condenses 
into a fog of which 95% of the particles are between 
0.2 jum and 0.22 ^t/m in diameter. This mean seed parti- 


cle size is small enough to minimize particle tracking er- 
rors for the vortex strengths found in these experiments 
(Ref. 33). The fog/air mixture was passed through a 
series of ventilation ducts and introduced into the rotor 
flow field at various locations. The entire test area was 
uniformly seeded by running the fog generator before 
each measurement. The 30,000 ft 3 (850 m 3 ) per minute 
capacity of the seeder allowed the entire test cell to be 
uniformly seeded in approximately 30 seconds. 

Signal bursts from seed particles passing through the 
measurement volume were received by the optics, and 
transmitted to a set of photomultiplier tubes where they 
were converted to analog signals. This analog signal 
was bandpass filtered from lOMHz-lOOMHz to remove 
the signal pedestal and any high frequency noise. The 
large range of the bandpass filter was required to allow 
measurement of the large flow reversal associated with 
the convection of a vortex core and the blade passage. 
The analog signal was digitized and sampled using a 
digital burst correlator. The flow velocities were then 
converted into three orthogonal components based on 
measurements of the beam crossing angles. Each mea- 
surement was phase-resolved with respect to the rotating 
blade by using a rotary encoder, which tagged each data 
point with a time stamp. The temporal phase-resolution 
of the encoder was 0.1°; however, the measurements 
were averaged into one-degree bins. The measurement 
uncertainty in using a one-degreebin width has been dis- 
cussed by Martin et al. (Ref. 24). 

Results and Discussion: 

Tip Shape Effects 

Phase-resolved LDV measurements documenting the 
formation, roll-up, and evolution of the blade tip vor- 
tex were acquired by making a radial traverse across the 
vortex core at various planes below the rotor. This al- 
lowed the core properties to be studied as a function of 
wake age, (see Fig. 4). The measurements for three 
other tip shapes were compared with the baseline blade. 
The planforms included a 20° swept back tip, a tapered 
tip with a 0.4 taper ratio, and a subwing tip. The vor- 
tex properties estimated from the velocity profiles were 
compared to the baseline measurements. In addition to 
the comparison, each of these test cases can be used for 
validating the ability of aerodynamic models to repre- 
sent the effects of tip shape on the wake geometry and 
vortex properties. 

Baseline Blade Vortex Properties 

The vortex measurements for the baseline blade are re- 
produced for comparison from Ref. 23. Figure 4 shows 
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Figure 3: Rotor test stand and LDV system. 

the measured wake geomtery. Figure 5 shows the varia- 
tion of the vortex core radius with wake age, and Fig. 6 
shows the variation of peak swirl velocity. In both of 
these figures, the aperiodicity corrections (based on light 
sheet measurements of core position) are also compared 
to the uncorrected measurements. 



Figure 4: Schematic showing measured wake geometry. 



Figure 5: Vortex core radius (inferred from velocity pro- 
files) as a function of wake age for the baseline blade. 


Wake Geometry 

The test case parameters for the swept and tapered tip 
blades are shown in Tables A2 and A3, respectively. In 
both cases, the blade was of constant chord from the 
root cut-out at 20% radius to the break-point at 80% 
radius. Estimates of the bound circulation distribution 
are shown in Fig. 7(a) using blade element momen- 
tum theory (BEMT) (Ref. 41). For comparison, the 
phase-averaged inflow measurements are shown across 
the TPP in Fig. 7(b). Note that without a force balance, 
the peak bound circulation (and therefore the tip vortex 


strength) could not be matched exactly. Even so, each of 
these test cases provides a useful set of measurements to 
compare with numerical predictions. 

One of the primary differences shown by the change 
in tip shape was the measured wake geometry shown in 
Fig. 8. The primary effect of taper is to decrease the 
bound circulation near the tip [see Fig. 7(a)], and this 
also decreases the effective downwash behind the blade 
[see Fig. 7(b)], As a result, the vortex tends to stay near 
the TPP for a larger fraction of a rotor revolution while 
moving radially inward at a much faster rate, as shown 
in Fig. 8. The primary effect of sweep appears to be that 
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Figure 6: Vortex peak swirl velocity as a function of 
wake age for the baseline blade. 




Figure 7: (a) Bound circulation distribution predicted us- 
ing BEMT, (b) Measured phase-average inflow distribu- 
tion. 

the vortex core is trailed outboard of the 100% radial sta- 
tion. The reason for this is shown by a top view in Fig. 9, 
where the outer edge of the swept tip does not conform 
to the circular streamlines. This tip shape is actually a 
“sheared” tip, with the airfoil chordlines perpendicular 
to the unswept line of aerodynamic centers (as opposed 
to the swept line of aerodynamic centers for true sweep). 
The final result is that sweep appears to decrease both 
the radial and axial convection of the vortex core, while 


Figure 8: Measured wake geometry for different tip 
shapes. 



Figure 9: Effect of swept tip in moving the initial core 
outside of r jR — 1 . 


taper appears to increase the radial convection and de- 
crease the axial convection. 

Vortex Core Properties 

The comparison of different tip shapes using vortex ve- 
locity field measurements requires great care in the ex- 
perimental methods. Preliminary measurements were 
used to identify the locations of the vortices and to un- 
derstand issues related to measurement uncertainty and 
spatial resolution. The variation of the vortex core radius 
(estimated from the velocity measurements) is shown in 
Figs. 1 0 and 1 1 as a function of wake age for the prelim- 
inary measurements and also for the refined measure- 
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ment technique, respectively. Over the first quarter of 
a rotor revolution, the vortex core size is on the order 
of the boundary layer thickness leaving the trailing edge 
for each of the blades. The error in the preliminary mea- 
surements over this region is significant when compared 
to the refined technique. In addition, clear trends and 
vortex structure are difficult to define in the presence of 
large measurement scatter. It is also important to note 
that as the spatial resolution degrades at early wake age, 
the variation in core radius with wake age appears con- 
stant. 

There is little difference in the vortex core growth 
shown in Figure 10 between the various tip shapes. As a 
result, the baseline aperiodicity measurements in Fig. 5 
are assumed to be a good approximation for the other tip 
shapes - with the exception of the subwing tip. In terms 
of the peak swirl velocity, differences appear for the var- 
ious tip shapes. A comparison is shown in Fig. 12, where 
the effects of blade loading variation between the cases 
have been removed to a first order by normalizing the 
peak swirl velocity by the average lift coefficient. The 
reason for this normalization is that the main focus of 
tip shape modification for B VI reduction is to minimize 
the peak swirl velocity for a given blade loading. 

The vortex peak swirl velocity (divided by tip speed)is 
a function of three factors, and these are the nondimen- 
sional vortex strength, the ratio of core circulation to the 
far-held value, and the nondimensional vortex core ra- 
dius. 



The vortex strength is directly proportional to the peak 
bound circulation, and this is the strongest influence on 
peak swirl velocity. 



where k\ is the ratio of the vortex strength to the peak 
bound circulation 

*■-£ ® 
Typical values appear to be k\ = 1 for an untwisted 
blade, and k\ — 0.7 for a highly twisted (tiltrotor) blade, 
but this needs additional research. The peak bound cir- 
culation is a function of blade planform (twist, taper, and 
sweep) and blade loading. 



where the proportionality constant, ki, is defined by 
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Figure 10: Variation of vortex core radius (uncorrected) 
with wake-age (a) Preliminary measurements (b) Mea- 
surements using improved technique. 


Typical values for ki in hover are 2 for an ideally twisted 
rectangular blade and 3 for an untwisted rectangular 
blade. 

The primary influence of tip shape is found in the last 
three terms of Eq. 4. By diving the peak swirl velocity 
by the factor 6Cr/c (average lift coefficient), the first 
term of Eq. 4 is removed. This allows the comparison of 
the tip shapes to avoid the influence of slightly different 
values of blade loading. This is shown by comparing the 
results in Figs. 1 1 and 12. 

A primary example of this reasoning is shown by the 
tapered tip blade. The estimated peak bound circula- 
tion for both the swept and tapered tips is approximately 
equal (see Fig. 7); however, the blade loading of the ta- 
pered tip is slightly higher. As a result, the tapered tip 
blade reduces the initial peak swirl velocity by reduc- 
ing the peak bound circulation (i.e., reducing k?) for a 
given blade loading. This benefit appears to decrease 
with wake age. It is interesting to note that once the dif- 
ference in blade loading is accounted for, the swept and 
baseline blades show similar behavior. 
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Figure 1 1 : Variation of vortex peak swirl velocity (un- 
corrected) with wake-age (a) Preliminary measurements 
(b) Measurements using improved technique. 
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Figure 12: Variation of vortex peak swirl velocity nor- 
malized by blade loading with wake-age. 


Subwing Blade 

The subwing blade concept has been investigated in 
wind tunnels (such as part of the NASA XV- 15 test pro- 
gram in Ref. 30) as well as in flight tests. A detailed 
survey of sub wing rotor blade tests is given by Ref. 25. 
In this current test, the subwing blade was constructed 
starting with a normal rectangular blade. Based on the 
results of Tangier (Ref. 31), the subwing chord was de- 
signed to be about 20% of the main blade chord and was 
placed at the quarter-chord position on the blade. The 
profile of the subwing was a fiat plate set at a pitch an- 
gle the same as that of the main blade. Manufacturing 
constraints resulted in a subwing tip having a chord of 
19.3% of the main blade chord, with the quarter-chord 
of the subwing being aligned with the 26.4% chord posi- 
tion on the main blade. The operating parameters for the 
subwing experiments are shown in Table A4, and were 
meant to be nominally the same as those used for the 
baseline experiments. As the subwing added less than 
3 % to the total blade area, its addition was assumed to 
cause no difference in the geometric or thrust- weighted 
rotor solidities. A single teetering blade was used and 
was balanced by a counterweight. 

As with the other blades, the phase-averaged inflow 
was used to provide an estimate of the thrust coefficient 
of the rotor. The inflow was measured across the TPP, 
as shown in Fig. 13. Based on these results, BEMT 
was used numerically to find a rotor thrust coefficient of 
0.0033. Using the blade solidity, the blade loading coef- 
ficient was then determined to be 0.0959. This is signifi- 
cantly higher than that found for the baseline blade, and 
was unexpected because of the nearly identical operating 
conditions (blade pitch and RPM). The exact cause of 
this dissimilarity is not currently understood; however, 
changes in performance caused by the subwing design 
have been documented by several other tests (these are 
summarized in Ref. 25). One possible reason is that the 
subwing reduces the tip loss of the main blade as evi- 
denced by the difference in the phase-averaged inflow 
near the tip. 


Subwing Blade: Flow Visualization 

Visualization of multiple vortex structures from unique 
tip shapes such as the subwing tip (Ref. 25) improve the 
understanding of how the tip shape can be used to mod- 
ify the resulting tip vortex structure. For comparison, the 
flow visualization images were taken at the same wake 
ages as were used for the baseline blade. The result- 
ing images are shown in Figs. 15 through 21, covering 
a wake age range from 10° to 209°. As with the base- 
line blade flow visualization, the scale is the same in the 
all the images. Figure 14 is a schematic of the first flow 
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Figure 13: Phase-averaged inflow velocity measured in 
tip-path-plane of the subwing blade compared to the 
baseline. 

r/fi | 



Figure 15: Flow visualization of wake generated by the 
subwing tip blade at £ = 10°. 



0.5 chord 



Figure 14: Schematic of the flow visualization image 
area for the sub wing blade (only the outer 12% of main 
blade radius is shown, and the direction of rotation is out 
of the page). 



visualization image at £ = 10°, and gives an idea of the 
scale of the images. It can be seen in Fig. 14 that the 
sub wing extends almost 4% radius outboard of the nom- 
inal blade radius, and the visible part of the blade corre- 
sponds to the outer 12% of the blade tip. 

The first picture (Fig. 15) shows the sub wing vortex 
immediately after formation, as the blade is rotating out 
of the plane of the page. The subwing vortex appears to 
originate at about the 80% radial station on the subwing. 
This was confirmed by slowly varying the phase of the 
image during the experiment. At this early wake age, 
the primary vortex from the main blade is not visible 
because optical access is blocked by the blade tip. 

Figure 16 shows the trailed vortices at £ = 25°, where 
the extreme tip of the blade can be seen on the far right. 
The vortex pair has rotated clockwise, with the subwing 
vortex moving axially up and radially inboard, and the 
primary vortex moving axially downward and slightly 
radially inboard. There is a clear difference in the ap- 
pearance of the vortices, with the subwing vortex show- 
ing a well-defined seed void, while the primary vortex 
shows a less clearly defined region at its center. This in- 


Figure 16: Flow visualization of the sub wing tip blade 
at £ .== 25°. 

dicates a relatively more diffuse primary vortex from the 
main blade compared to that produced by the subwing. 

The vortices at a wake age of 45° are shown in Fig. 17. 
At this point, the pair has rotated in the light sheet plane 
about 100° from the original position, and the subwing 
vortex is now above the primary vortex. The vortex sheet 
can be seen to wrap up around the outboard side of the 
primary vortex. 

The image in Fig. 18 shows the vortices at £ = 67°. 
Here the vortex pair has rotated about 125° in the light 
sheet plane, and the primary vortex is now outboard of 
the subwing vortex. The subwing vortex lies in the tip- 
path-plane. At this point the two vortices can be seen 
to be beginning their merging process, as their struc- 
tures seem to be undergoing some kind of interaction. 
The aperiodicity of the locations of the vortex cores was 
noted to increase at this point. 
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Figure 18: Flow visualization of the sub wing tip blade 
at C = 67°. 


Figure 19 shows the vortices at £ = 90°, where they 
have rotated about 180° in the light sheet plane and have 
now switched positions relative to their initial orienta- 
tion. The seed voids are present, but there is obviously 
some interaction occurring between the vortices. The 
primary vortex, in particular, appears very diffuse and 
almost no clear structure can be seen. The vortex sheet 
is seen to still be in the process of being entrained into 
the primary vortex. The turbulent flow from the vortex 
sheet seems to be creating significant instability within 
the primary vortex and in the surrounding area. 

At a wake age of 137°, as shown in Fig. 20, the vor- 
tices appear to have completed merging and are now es- 
sentially one tip vortex. However, the aperiodicity of the 
vortex core location was noted to still be significant at 
this wake age. A seed void is visible but it seems very 
diffuse. In fact, the entire vortex appears to be more 


Figure 20: Flow visualization of the subwing tip blade 
at £ = 137°. 


diffuse, and is surrounded by flow that appears to con- 
tain turbulent eddies. The vortex sheet has convected 
downward faster than the vortex core, and only the out- 
board edge of the sheet can be observed. However, the 
sheet appears to still be wrapping into the merged vor- 
tex, causing more turbulent flow to be entrained into the 
flow surrounding the core. 

Results for the last wake age of 209° are shown in 
Fig. 21. No clearly laminar flow region can be seen. 
The vortex sheet has wrapped all the way around into the 
merged vortex. Note that as mentioned previously, these 
images are all to the same scale. The very large size of 
the seed void and of the overall vortex structure at this 
wake age confirms the increase in core radius expected 
from the merged vortices generated by the subwing tip 
blade. 
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Figure 21: Flow visualization of the sub wing tip blade 
at 5 = 209°. 

Subwing Blade: Measurement Technique 

The creation of multiple vortices required a different 
technique to capture the vortex velocity profiles. Using 
the flow visualization results for the subwing blade, it 
was possible to construct radial measurement grids for 
the phase-resolved LDV measurements. These grids are 
shown in Fig. 22, along with a scale representation of the 
the subsequently estimated vortex core sizes as detailed 
in Ref. 26. As with the baseline blade, the measurement 
grids were designed to extend between 5 and 10 core 
radii on either side of the vortices. Because of the ex- 
tremely small physical size of the subwing vortices, the 
grids were constructed with center portions having spac- 
ing of 0.25 mm (0.5 %c) between measurement points so 
as to accurately resolve the large velocity gradients sur- 
rounding the vortex core (the variable grid spacing is not 
shown in Fig. 22). 

At the earliest wake ages, the difficulty of finding the 
exact location of the very small subwing vortex necessi- 
tated using a vertical measurement grid (parallel to the 
rotor shaft), rather than the standard radial grid. This 
ensured that the grid would capture the subwing vortex 
regardless of where it crossed vertically. Note that be- 
cause two of the measured cases used vertical grids, their 
results cannot be plotted on the same axes as the other 
wake ages, and thus they are shown separately. How- 
ever, to enable comparison, the plots are all shown with 
the same scale and orientation. This allows the results 
for the tangential and axial velocity components to be 
directly compared. 


Subwing Blade: Vortex Trajectories 

A closer view of the wake geometry in terms of the axial 
or vertical (z) and radial (y) locations of the vortices at 



Figure 22: LDV measurement grids for the subwing 
blade. 

each wake age measured is shown in Fig. 23. The wake 
ages measured with the subwing tip were from £ — 8° to 
£ = 260°. It was found through experimentation that be- 
tween about £ = 40° and £ = 130° the interaction and 
merging of the subwing and primary vortices created 
a less coherent vortical structure that could not be ac- 
curately measured because of excessive aperiodicity in 
the flow. While the flow at these wake ages could still 
be qualitatively analyzed through flow visualization, the 
statistic ally-based measurement nature of LDV made it 
extremely difficult to obtain meaningful results in this 
region. After the merging region, the vortex could again 
be measured with aperiodicity levels of the same order as 
the baseline. Detailed aperiodicity measurements were 
not made, however, as in the case of the baseline blade. 

It can be seen that the subwing vortices at the two ear- 
liest wake ages (^ = 8° and 22°) lie outboard of the blade 
tip (where y/R = 1 is the tip of the main blade), and at 
the second wake age (£ = 22°) the vortex additionally 
lies above the tip path plane (i.e., z/R < 0). This is a 
result of the rotation of the subwing vortex around the 
stronger main tip vortex. 

Subwing Blade: Tangential Velocity Profiles 

The vortex-induced tangential velocity profiles, normal- 
ized by the rotor tip speed, are shown in Figs. 24 and 
25. The subwing vortex and the primary vortex were 
studied independently, because only one vortex was tar- 
geted per LDV traverse across the flow field. One of 
the most significant contributions of this work is simply 
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Figure 23: Sub wing blade vortex locations relative to 
tip-path-plane. 

the process of capturing these highly dynamic interact- 
ing vortices before and after the merging process. It is 
important to note that before the merging process there 
are two vortices present in the flow, so neither vortex 
can be isolated from the induced effects of the other in 
these types of experiments. For this reason, asymme- 
tries in the profiles are very different than in the case 
of the baseline blades. It is not very accurate to deduce 
from the measurements any vortex properties (e.g., core 
radius and peak swirl velocity) that inherently assume 
a single isolated vortex. In fact, it is the interaction of 
the induced velocity fields that drives the merging pro- 
cess. Keeping this limitation in mind, results are shown 
in an effort to provide best estimates of these properties. 
The measurement grids shown in Fig. 22 were chosen, 
in part, to minimize the signature of the adjacent vortex. 
The purpose is to understand the dynamics and timing 
involved in the merging process. The most accurate re- 
sults (and the most practical) are the final vortex proper- 
ties and structure after the merging process. The reason 
is that the final merged vortex structure is what (in most 
cases) would be encountered by a following blade. 

Subwing Blade: Axial Velocity Profiles 

The axial velocity profiles, normalized by the tip speed, 
for each of the six measured wake ages are shown in 
Figs. 26 and 27. Again the subwing and primary vortices 
were measured with different grid orientations, and so 
are shown separately. 


Figure 24: Measured tangential velocity profiles of sub- 
wing vortices at two wake ages. 
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Figure 25: Measured tangential velocity profiles of pri- 
mary and merged vortices at four wake ages. 

It can be seen from Fig. 27 that at the early wake ages, 
the primary vortex has axial velocities of comparable 
magnitude to the tangential velocities. From the stabil- 
ity analysis of Ref. 26, the stability of the vortex core 
is a function of the ratio of peak swirl velocity to peak 
axial vleocity deficit. In this case, both cores would be 
considered unstable according the peak values shown in 






Figure 26: Measured axial velocity profiles of sub wing 
vortices at two wake ages. 



Figure 27: Measured axial velocity profiles of primary 
and merged vortices at four wake ages. 

Figs. 25,24,26, and 27. This is appears to be consistent 
with the bands of turbulent flow surrounding both cores 
observed in the flow visualization results. The axial ve- 
locity peak rapidly decreases after the merging process 
is complete. The width of the axial velocity deficit re- 
gion appears to be 2 to 3 times the core radius on each 


side of the vortex, as in the case of the baseline blade 
In terms of axial and tangential velocity profiles, the fi- 
nal structure of the core appears to be very similar to the 
baseline blade once the merging process is complete 

Subwing Blade: Vortex Properties 

The vortex core properties can be estimated from the ve- 
locity measurements, keeping in mind the above men- 
tioned limitations concerning the difference between an 
isolated vortex and a pair of vortices. One of the primary 
effects of the subwing blade is a significant increase in 
the net core size after the merging process. Figure 28 
shows the change in estimated viscous core radius with 
wake age, normalized by the blade chord. As with the 
baseline blade, the vortex core radius is defined to be 
half the distance between the measured peaks in the tan- 
gential velocity. 

At the earliest wake age of £ « 10°, the core radius 
of the subwing vortex was estimated to be 2.5% of the 
blade chord, while the core of the primary vortex was 
7.7% chord. The core grows to about 23% chord after 
the vortices have merged, and stays nominally constant 
thereafter. By comparing the subwing blade with the 
baseline, the core radius appears to be approximately 
double just after the merging process is complete. A 
large increase in core radius would be an advantage for 
alleviation of noise and vibratory airloads, especially if 
the peak swirl velocity decreased by the same factor. 

These measurements show this might not be the case, 
however, because the merging process appears to redis- 
tribute the vortex circulation distribution. In the flow 
visualization and the measurements, both the subwing 
and primary vortex appear to rapidly diffuse before the 
merging. Accurate modelling of this type of blade tip 
requires the ability to capture both the individual diffu- 
sion of each vortex core, and the global merging pro- 
cess. The entire process appears to be the result of a 
viscous-inviscid interaction. Even these basic measure- 
ments pose a tough challenge for current modeling ca- 
pabilities. 

This redistribution of circulation is shown by the 
variation of non-dimensional peak tangential velocity, 
Vtomax/^R' with wake age, shown in Fig. 29. For 
comparison, the baseline measurements (uncorrected) 
are also shown in this figure. The primary difference 
between the baseline and the subwing blades is that the 
subwing blade has a very small change in peak swirl ve- 
locity with wake age. The reason for this is an inviscid 
-viscous interaction that results in almost no net change 
in the swirl velocity as the vortex ages. 

The interaction is described as follows. The total vor- 
tex strength trailed behind the blade is divided between 
the two vortices. When the vortices merge, the effective 
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Figure 28: Vortex core radius versus wake age for the Figure 29: Vortex peak tangential velocity versus wake 
subwing blade compared to the baseline. age for the subwing blade compared to the baseline. 


circulation approximately doubles from the value each 
core showed before merging. This is the inviscid part 
of the interaction - which can be predicted by poten- 
tial flow vortex conservation laws. The viscous part of 
the interaction is that the effective core size increases 
by roughly a factor of two from that of the largest core 
before merging. It was shown in Ref. 26 that the peak 
swirl velocity is by definition proportional to the ratio 
of the vortex circulation to core radius. A doubling of 
the circulation and a doubling of the effective core ra- 
dius cancel. This results in a very small change in the 
peak swirl velocity after merging. In these tests the cir- 
culation is not equally divided, and so there is a reduc- 
tion in the peak swirl velocity when the comparison is 
made by normalizing with blade loading (as shown by 
Fig. 30). The benefit in terms of reduced peak swirl ve- 
locity appears to decrease with wake age, however, more 
measurements are required. 

These results are supported by the flow visualization 
images at £ = 209°. The vortex from the baseline blade 
[Fig. 3 1(a)] has a clear seed void at its center, surrounded 
by a small region of relatively smooth stratified flow. 
Further out, the flow becomes much more turbulent, with 
eddies being visible. The merged vortex from the sub- 
wing tip blade has a seed void, but it is much less clearly 
defined [Fig. 31(b)]. Also, nowhere does the flow ap- 
pear to be smooth and laminar. Turbulent flow is seen to 
completely surround the vortex core. 



Figure 30: Vortex peak tangential velocity normalized 
by blade loading versus wake age for the subwing blade 
compared to the baseline. 

Subwing Blade: Scale Effects 

From the subwing results, the merging process is as- 
sumed to be a strong function of the vortex Reynolds 
number. For this reason, the same type of increase in 
core radius might not be observed at full-scale flight 
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(a) Baseline blade wake 



Figure 31: Comparison of flow visualization of wake at 
C = 209°. 


provided insight into the basic fluid mechanics. These 
results are presented here in terms of vortex structure 
and vortex core stability. An understanding of both of 
these topics supports the development and analysis of 
tip vortex alleviation devices. 


Iverson’s Equation 

Iverson (Ref. 34) developed a correlation function based 
on the self-similar turbulent decay of a line vortex. The 
objective of developing a correlation function was to un- 
derstand the scale effects involved in measuring vortices 
in low Reynolds number facilities compared to the ac- 
tual high Reynolds number behavior measured in flight 
tests. The purpose of solving Iverson’s equation is to de- 
termine a similarity solution for the vortex velocity pro- 
files for the case when a turbulent viscosity is allowed 
to vary across the vortex. The importance of obtaining 
a similarity solution is that it is an exact solution of the 
Navier-Stokes equations, which is useful for developing 
fundamental vortex models such as the Lamb vortex. A 
unique aspect of Iverson’s result is the Reynolds number 
dependence on the final solution, which does not occur 
in the case of the Lamb vortex. 

The vorticity transport equation with a variable eddy 
viscosity was given by Iverson in terms of circulation 


ar a 

a t ~ r dr 



a 

+ 2v r r— 
ar 



( 6 ) 


where the turbulent viscosity is a function of the dis- 
tance from the center of the core and the local circulation 
strength 


Vj = 


a 2 r 2 
2n 


r dr 


+V 


(7) 


Reynolds numbers. More research is needed to under- 
stand the vortex structure that would be formed in for- 
ward flight. These tests, however, show that it is an ad- 
vantage to divide the circulation between two vortices in 
terms of the magnitude of the initial peak swirl velocity. 
It appears the benefit is reduced as the vortices merge 
in comparison with the baseline blade. It is concluded 
in terms of BVI noise reduction, that designs aimed at 
dividing the circulation while preventing merging might 
be an advantage over the subwing blade - especially at 
late wake ages where the most severe parallel interac- 
tions occur. 


Analysis of Baseline Vortex 
Measurements 


This is equivalent to specifying a variable mixing length 
as 

L = ar (8) 

where the empirical constant according to Iverson is 

a = 0.01854 (9) 


This equation is non-dimensionalized in order to ob- 
tain a single equation that provides a similarity solution 


3r| 2 


2x 1 \ _ dy 

^3t| Y a 2 Re v ) 3t| 


with the similarity variables 




( 10 ) 


( 11 ) 


The baseline blade measurements were first reported in This equation will be referred to as Iverson s equation. 

Ref. 23. Recently, more detailed analysis of the data has The solution (outlined in Ref. 26) provides a self- similar 


2072 




Figure 32: Vortex core strength as a function of 

Reynolds number obtained by solving Iverson’s equa- 
tion compared to several sets of measurements compiled 
in Ref. 26. 



Effective lift coefficient, C L or 6C T /a 

Figure 33: Vortex initial peak swirl velocity as a func- 
tion of effective lift coefficient for rotor-blades and non- 
rotating wings (from Ref. 26). 


circulation profile, y(r|) across the vortex core. This type 
of profile helps to describe the structure of the vortex. 

The core strength y(l) = r c /r v is found to be a func- 
tion of vortex Reynolds number, as shown in Fig. 32. 
At low Reynolds numbers the fully laminar solution is 


identical to the Lamb vortex. At high Reynolds numbers 
the fully turbulent solution is another fundamental exact 
solution for a vortex flow. As shown by the results in 
Fig. 32, most experiments have been made in the transi- 
tional range of the solution. The advantage of Iverson’s 
method is that it provides some fundamental self-similar 
profiles to use for comparing experimental data. The dis- 
advantage is that the prediction of the core strength does 
not appear to agree with available measurements. This is 
shown by the comparison between the theory and mea- 
surements in Fig. 32 of the ratio r c /F v . One possible 
reason for this is that 2-D vortex models based on a so- 
lution of the Navier-Stokes equation do not take into ac- 
count the variation of the initial core structure with lift 
coefficient (i.e., the initial conditions at fo vary with the 
blade lift). This is discussed by McCormick (Ref. 35), 
and is supported by the variation of initial peak swirl ve- 
locity with lift coefficient shown in Fig. 33. 

Furthermore, other factors may influence the role of 
turbulence in the core structure, such as tip separation 
effects, boundary layer states near the trailing edge, and 
roll-up of the trailing vortex sheet. For this reason, air- 
craft vortices might not be properly modeled by the ideal 
laminar or ideal turbulent solutions for an isolated two- 
dimensional vortex. 


Extension of Iverson’s Result to Transitional 
Vortices 


The baseline measurements in this study appear partly 
laminar and partly turbulent (see Ref. 26), and so these 
two fundamental solutions are superimposed to obtain a 
model for transitional vortices. The important question 
is where the vortex transitions from laminar to turbulent, 
and the ability to predict this behavior. A preliminary 
answer to this question can be found using the concept 
of relaminarization and the Richardson number formu- 
lation presented in Ref. 26. 

Using the definition of a shape factor applied to Iver- 
son’s equation 

5 - .. YllL ~ 1 271 ( 12 ) 

dVg/dr CnAlc) yin (\AlAlc) + 2dy/dr\ 
the Richardson number is 


Ri = 2S p (S p + l) (13) 


The procedure begins with obtaining three solutions for 
the vortex using Iverson’s method. One solution is fully 
turbulent y r, one is fully laminar y t, and one is at the 
given vortex Reynolds number y (Re v ). 

The next step is to test for relaminarization, where the 
test parameter is defined as (Ref. 26) 


Rl = 


Ri 

Rel /4 


(14) 






Figure 34: Distribution of Richardson number across 
vortex in terms of Rl . 

When Ri < 1 no relaminarization can occur, and when 
Ri > 1 the streamline curvature limits the role of turbu- 
lence in transporting vorticity. 

Figures 34 and 35 show the distribution of Rl outside 
of the core for the three cases. In the case of the fully tur- 
bulent vortex, Rl remains less than one, and in the case 
of the fully laminar vortex, Rl remains greater than one 
over the region where 99% of the vorticity is contained. 
For the case of Re v = 50800 (which corresponds to the 
baseline test case), the value of Rl becomes less than 1 
near r/r c = 2. Using the value of at the core bound- 
ary, a transitional circulation profile y rr is constructed by 
superposition of the fully laminar y l and fully turbulent 
Jt profiles. The superposition is accomplished using a 
form of the Gudermannian function 

( ‘ 5) 

The transitional circulation profile is then 

Yfr=/Yr + (1 -/)Yi 06) 

and the corresponding velocity profile is shown in com- 
parison with the measurements in Fig. 36. The result is 
in agreement with the measurements. 

Repeating this procedure over a range of Reynolds 
numbers creates a family of transitional vortex profiles. 
The only unknown is the radial position where 50% of 
the transition takes place, \Rl c - The variation of Rl c 
with Reynolds number is shown in Fig. 37. This figure 
shows the range of Reynolds numbers where the vortex 
might appear transitional. More research is needed in 
this area. 



Figure 35: Transition point defined by relaminarization 
test. 
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Figure 36: Transitional vortex velocity profile Re v = 
50800 compared to laminar and turbulent profiles and 
measurements, Re v — 50800. 


Vortex Core Stability 

Another decay mechanism that may also participate in 
the core growth of rotor vortices is the stability of the 
vortex core. This is related to axial flow and can be 
described by an adaptation of the Leibovich-Stewartson 
vortex stability theory (Ref. 36, 37). This practical the- 
ory is applied by first fitting the following Lamb-like 
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Figure 37: The variation of Ri c with Reynolds number 
obtained from the solution of Iverson’s equation. 


functions to the measured velocity distributions using 


V%_ 

SIR 


Vox 

SIR 



(17) 

Ce~ D ' 2 

(18) 


where the unknown coefficients are A, B, C, and D. The 
stability parameter is then defined by 


SIR \BCj 


(19) 


and the vortex core is considered unstable when S < 0.9 
(Ref. 38). Using this criteria, the present measurements 
of peak axial and swirl velocity show that S = 0.65, so 
that initially the vortex at £ = 3° may be unstable, and 
by a significant margin. This is a sufficient condition 
for the vortex to undergo a “hydraulic jump” to a larger 
diameter, or even burst (Ref. 38). The “hydraulic jump” 
process is described by Rule (Ref. 39), and is supported 
by the measurements of the changing core circulation 
structure in Ref. 26. The strength of this instability is 
proportional to the ratio of the peak swirl velocity to the 
axial velocity deficit. 

The Leibovich-Stewartson vortex stability theory con- 
siders only the general structure of the entire core. The 
final conclusion is that the core is either stable or unsta- 
ble, but there is no indication of how the different re- 
gions of the vortex core contribute to the stability. De- 
tailed modeling of the vortex, as well as effective allevi- 
ation strategies, requires an understanding of which re- 
gions of the vortex are the most stable and unstable. This 


Figure 38: Self-similar tangential velocity distribution 
compared to transitional Iverson model. 



Figure 39: Self-similar axial velocity distribution esti- 
mated from early wake age measurement at £ = 3°. 

is very similar to the analysis of boundary layers, where 
the mean flow properties are used to determine a shape 
factor, H, that allows the regions of laminar and turbu- 
lent flow to be estimated. For vortex flows the same type 
of procedure involves using the Ludwieg stability test 
(Ref. 40) to classify regions of the core as either stable 
or unstable in terms of the ability to produce or damp out 
turbulence. 
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Figure 40: Ratio of peak axial velocity to peak swirl 
velocity as a function of wake age. 



Figure 41: Ludweig’s stability coefficients as a function 
of the radial position across the vortex. 


The procedure for calculating the Ludwieg stability 
map of the vortex core requires both the tangential and 
axial velocity profiles. These have been measured in the 
current experiments and can be used along with the the- 
ory to estimate the vortex structure from the mean flow 
properties. The Ludwieg stability theory involves two 
coefficients that are simply a measure of the ratio of the 
local velocity gradient to the local angular velocity of the 
fluid. One coefficient is determined from the tangential 
velocity distribution 




( 20 ) 


where all of the variables are assumed to the self-similar 


values determined from the measured profile shown ' 
Fig. 38. The other coefficient is determined from the 
axial velocity deficit measured at the earliest wake ao 
- see Fig. 39 



, __ (^c)q / dVaxjdr \ 

Z ( V e) ma x V V$/r ) ( 2l ) 


Note the sign convention of is that a deficit (wake- 
like profile) is a positive increment from the free-stream 
velocity. The maximum value at the center of the profile 
is (V^) 0 . The coefficient C z is also a function of the ratio 
of the peak axial velocity to the peak swirl velocity. 

The variation of the tangential and axial coefficients is 
shown in Fig. 41. In this figure, it was assumed that the 
ratio of the peak axial velocity to the peak swirl velocity 
was one, so that the results can be scaled for different 
wake ages by multiplying C z by the actual ratio shown 
in Fig. 40. In addition, a curve fit to both profiles was 
used to calculate the velocity gradients. The curve fit to 
the axial profile is shown as a simple exponential distri- 
bution in Fig. 39. The curve fit to the tangential velocity 
profile shown in Fig. 38 is the transitional Iverson result 
as derived. 

The final step after the coefficients have been deter- 
mined is to make a comparison with the stability bound- 
ary. The stability boundary is a function of these two 
coefficients, and Ludwieg’s condition for stability is de- 
fined (Ref. 40) by the inequality 



1 

i 


(l-C 6 )(l-Ci)-C?(S-C e 


> 0 ( 22 ) 


By plotting the values of C e versus C z for different radial 
positions from the vortex center, the stability of different 
regions of the vortex core can be estimated with respect 
to the stability boundary shown in Fig. 42. As shown in 
this figure, there is a band of radial positions where the 
vortex structure is estimated to be unstable, and there- 
fore capable of producing turbulence. This appears to 
correlate well with the flow visualization image shown 
in Fig. 44. As the ratio of the peak axial velocity to peak 
swirl velocity decreases with wake age (Fig. 40), the 
measured curve would scale below the stability bound- 
ary and the vortex would become stable. 

A more physical presentation of the same result can 
be obtained by plotting the difference between the mea- 
sured results and the stability boundary, as shown in 
Fig. 43. The fundamental result obtained by applying 
the Ludwieg stability test to measurements at the ear- 
liest wake age of 3° is that the vortex core is stabi- 
lized by the rotational flow. These results show that if 
there is a turbulence producing region of the vortex, it 
would appear outside the core boundary from rf r c = 1 .5 
to r/r c = 2.75. This result is important because most 
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Figure 42: Ludweig’s stability map. 



Figure 43: Ludweig’s stability test as a function of posi- 
tion across the vortex core. 


ff vortex alleviation devices attempt to promote turbulence 

i inside the vortex core where instabilities are quickly 

V damped by the high streamline curvature. By using de- 

tailed measurements and stability analysis, future vortex 
alleviation devices might be designed to take advantage 
v of the unstable region outside of the core. In theory, such 

a device would not have the severe drag penalty of typi- 
cal devices designed to increase the axial velocity deficit 
inside the core. More research is needed to confirm this 

S V hypothesis. 



Figure 44: Flow visualization image zoomed-in on the 
vortex core region. 


Conclusions 

Three-dimensional phase-resolved laser Doppler ve- 
locimetry measurements with high spatial and temporal 
resolution were made in the viscous wake trailing from 
a helicopter rotor blade. Detailed velocity profiles of the 
trailing vortex filament enabled the study of the vortex 
core formation, vortex roll-up, and vortex evolution. The 
conclusions drawn from this study are as follows: 

1. One of the primary differences shown by the 
change in tip shape was the measured wake geom- 
etry. Tip sweep appeared to decrease both the ra- 
dial and axial convection of the core, while taper 
appeared to increase the radial convection and de- 
crease the axial convection. These results are for 
hover, and more research is needed to conclude the 
same effects would apply to forward flight. 

2. While there was little difference in the measured 
vortex core growth between the various tip shapes, 
differences were observed in the peak swirl veloc- 
ity. After the blade loading variation between tests 
is used to normalize the peak swirl velocity, the 
swept and baseline blades showed similar behav- 
ior. The effect of blade taper was to reduce the 
initial peak swirl velocity by a significant fraction. 
It appears that this is accomplished by decreasing 
the peak bound circulation (and therefore the trailed 
vortex strength) for a given blade loading. 

3. The sub wing measurements showed that between 
the wake ages of £ = 40° and C, = 130°, the in- 
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teraction and merging of the subwing and primary 
vortices created a less coherent vortical structure. 
In terms of axial and tangential velocity profiles, 
the final structure of the core generated by the sub- 
wing blade appears to be very similar to the base- 
line blade after the merging process is complete. 
The primary difference between the baseline and 
the subwing blades is that the subwing blade has a 
very small change in peak swirl velocity with wake 
age and a much larger core size. 

4. One source of vortex core instability is the ratio of 
the peak swirl velocity to the axial velocity deficit. 
These results show that if there is a turbulence pro- 
ducing region of the vortex structure, it is predicted 
to be outside of the core boundary from r/r c = 1.5 
to r/r c = 2.75. This is the region where the mea- 
sured vortex self-similar velocity profiles appear to 
transition from a laminar to a more turbulent model. 
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Appendix: Blade Parameters and 
Test Conditions 


Radius, R 

406.0 mm 

Chord, c 

44.5 mm 

Root cutout, r 0 

0.2R 

Tip speed, OR 

89.28 m/s 

Airfoil section 

NACA 2415 

Collective pitch, 0 q 

4.5° 

Twist, 0 rvv 

0° 

Rotational frequency, ft 

30Hz = 707t rad/s 

Tip Mach number, M t i p 

0.26 

Tip Reynolds number, Re t i P 

272,000 

Thrust coefficient, Cj 

0.0022 

Solidity, a 

0.035 

Blade loading coefficient, Cf/a 

0.064 

Table 1: Parameters of the baseline blade 


Radius, R 

407.99 mm 


Chord, c 

44.0 mm 

Root cutout, r 0 

0.2R 

Break point, r b 

0.8R 

Sweep, A 

20° 

Tip speed, OR 

89.72 m/s 

Airfoil section 

NACA 2415 

Collective pitch, 0 q 

4.5° 

Twist, Q tw 

0° 

Rotational frequency, ft 

219.9 rad/s 

Tip Mach number, M t f P 

0.3 

Tip Reynolds number, Rea p 

250,000 

Thrust coefficient, Cj 

0.00258 

Solidity, c 

0.034 

Blade loading coefficient, Ct/o 

0.0758 


Table 2: Parameters of the swept tip blade 


Radius, R 

410.00 mm 

Root chord, c 0 

44.0 mm 

Tip chord, c t 

17.6 mm 

Effective chord, c e 

37.1 mm 

Root cutout, r 0 

0.2R 

Break point, r b 

0.8R 

Taper ratio, z 

0.4 

Tip speed, QR 

91.36 m/s 

Airfoil section 

NACA 2415 

Collective pitch, 0 q 

6.0° 

Twist, 0 rvv 

0° 

Rotational frequency, ft 

219.9 rad/s 

Tip Mach number, M t i p 

0.3 

Tip Reynolds number, Re cip 

100,000 

Thrust coefficient, Cj 

6.00259 

Solidity, c e 

0.029 

Blade loading coefficient, Cr/o e 

0.089 


Table 3: Parameters of the tapered tip blade 


Radius, R 

407.0 mm 

Chord, c 

44.0 mm 

Root cutout, r 0 

0.2R 

Tip speed, OR 

89.504 m/s ' 

Airfoil section 

NACA 2415 

Collective pitch, 0q 

4.5° ’ 

Twist, 0 fw 

0° 

Rotational frequency, ft 

219.9 rad/s 

Tip Mach number, M t \ p 

0.3 

Tip Reynolds number, Re t i P 

250,000 T 

Thrust coefficient, Cj 

0.0033 ■; 

Solidity, a 

0.0344 

Blade loading coefficient, Cj/c 

0.0959 


• -.41 My. V T 

Table 4: Parameters of the subwing tip blade 





